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Abstract: In this paper we focus on the problem of the degree sequence for the following random 
graph process. At any time-step t, one of the following three substeps is executed: with probability 
qi, a new vertex xt and m edges incident with xt are added; or, with probability a — ai, m edges 
are added; or finally, with probability 1 — a, m random edges are deleted. Note that in any case 
edges are added in the manner of preferential attachment, we prove that there exists a critical point 
a c satisfying: 1) if ct\ < a c , then the model has power law degree sequence; 2) if ai > a c , then the 
model has exponential degree sequence; and 3) if a.\ = a c , then the model has a degree sequence 
lying between the above two cases. 

1 Introduction and statement of the results 

Recently there has been much interest in studying large-scale real-world networks and attempting 
to model their properties. For a general introduction to this topic, readers can refer to Albert and 



AMS classification: 05C 07. 05C 80. 

Key words and phrases: degree sequence; power law; critical phenomenon; real-world networks. 

* Supported in part by the Natural Science Foundation of China 

^Supported in part by the Natural Science Foundation of China under grants 10671197 and 10721101 

* Partially supported by the Natural Science Foundation of China under grants 60674082, 70221001 and 70731003. 
§ Supported by the Natural Science Foundation of China and MicroSoft Research Asia under grant 60633010 



1 



Barabasi pQ, Aiello, Chung and Lu [3], Bollobas and Riordan [9J, Hayes [T7], Newman [22] and Watts 
[25] , Although the study of real- world networks as graphs can be traced back to long time ago such 
as the classical model proposed by Erdos and Renyi [2] and Grilbert [IB], recent influential activity 
perhaps started with the work of Watts and Strogatz about the 'small-world phenomenon' published 
in 1998 [27]. Another influential work may be due to the scale- free model proposed by Bollobas and 
Albert in 1999 [5]. Since then various forms of scale- free phenomenon have been widely revealed. 
In particular, power law degree distributions have been extensively investigated. Many new models 
have been introduced to circumvent the shortcomings of the classical models introduced by Erdos 
and Renyi [14] and Grilbert [16j . One class of these new models was aimed to explain the underlying 
causes for the emergence of power law degree distributions. This can be observed in 'LCD model' [TO] 
and its generalization due to Buckley and Osthus [8j, 'copying' models of Kumar et al. [19], the very 
general models defined by Copper and Frieze [12] and the other model with random deletions defined 
by Copper, Frieze and Vera [13] etc. 

For the real- world network of World Wide Web/Internet, experimental studies by Albert, Barabasi 
and Jeong [2J, Broder et al. [7] and Faloutsos, Faloutsos and Faloutsos [15] demonstrated that the 
proportion of vertices of a given degree follows an approximate inverse power law, i.e., the proportion 
of vertices of degree k is approximately Ck~ a for some constants C and a. However other forms of 
the degree distributions can also be observed in real- world networks (see [3] and [25])- For example, 
Guassian distributions can be observed in the acquaintance network of Mormons [6j; exponential 
distribution can be observed in the powergrid of southern California [27]. On the other hand, the 
degree distribution of the network of world airports [4] interpolates between Gaussian and exponential 
distributions, whereas the degree distribution of the citation network in high energy physics [20] 
interpolates between exponential and power law distributions. For more forms of degree distributions, 
readers can refer to [24] . 

Different models often lead to different forms of degree distributions. An interesting problem 
arises naturally: does it exist some dynamically evolving random graph process which brings forth 
various degree distributions by continuous changing of its parameters only? This phenomenon has 
been numerically investigated in reference 28J: For a general model of collaboration networks in [28] 
, Zhou et al. indicate that, while a relevant parameter a increases from to 1.5, four kinds of degree 
distributions appear as exponential, arsy-varsy, semi-power law and power law in turn. Note that the 
above classification is rather rough as no unambiguous borderline between two neighboring patterns is 
determined. However, to the best of our knowledge, it seems that the problem and its answer have not 
been formulated in a mathematically rigorous manner. In this paper we focus on a model with edge 
deletions and provide precise analysis, while a parameter varies, the model exhibits various degree 
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distributions. 

Now, we begin to introduce our model and then state our main results. Consider the following 
process which generates a sequence of graphs G t — (V t ,E t ), t > 1. Write v t — \ V t \ and e t = \E t \. 

Time-Step 1. Let G\ consist of an isolated vertex x\. 
Time-Step t > 2. 

1, With probability a\ > we add a vertex xt to Gt-i- We then add m random edges incident 
with xt- In the case of et-\ > 0, the m random neighbours w\,W2, ■ ■ ■ , w m are chosen independently. 
For 1 < i < m and w € Vt—i, 

2e t _i 

where d w (t — 1) denotes the degree of vertex w at the beginning of substep t. Thus neighbours are 
chosen by preferential attachment. In case of &t-x — 0, then we add a new vertex xt and join it to a 
randomly chosen vertex in Vt_i. 

2, With probability a — a.\ > we add m random edges to existing vertices. If &t-\ > 0, then 
both endpoints are chosen independently with the same probabilities as in Otherwise, we do 
nothing. 

3, With probability 1 — a > we delete min{m, et-i} randomly chosen edges from Et—i. 

Remark 1.1 The deference between our model and the model introduced in is that, in our setting, 
vertex deletions, loop and multi-edge erasures are forbidden, which makes {et : t > 1} Markovian and 
makes it possible for us to give exact estimation to e t . 

In order to make the problem meaningful, the following inequalities are natural and necessary: 

1/2 < a < 1; < ai < a. (1.2) 

For given a and u\ satisfying (|1.2[) . define 

a c := Aa — 2, r\ := a c m/2, (1-3) 

and choose e = e(a, a±) G (0, 77) such that 

f m(a r — at ) 1 1 . 

p.i^max 1* <1. (1.4) 



2(7? -e) '2 



Note that in case of a\ > a c , p e = h. Let 



a a c a\-a c 2a c - a\ a c 

a c — ai 2(1 — a) 2a c 2(1 — a) 
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Obviously, j3 is well denned when a% ^ a c and < 7 < 1 when ot\ > a c . To get our main results, 
besides (|1.2[) . the following condition is necessary 



ol\ < 2a c 



(1.6) 



Now, Let Dk(t) be the number of vertices with degree k > in Gt and let Dk(t) be the expectation 
of Dk(t). The main results of this paper follow as 

Theorem 1.1 Assume that hl.2\) and (1.6\) hold. Then a c defined in hl.3\) is a critical point for the 
degree sequence of the model satisfying: 

1) if ai < a c , then there exists a constant C\ = Ci(m, a, a±) such that, for any v £ (0, 1 — p e ), 

D k (t) 



t 



-dk- 1 -? 



= 0(t p ' +v - l ) + 0{k- 2 -P)- 



2) if (Xi > a c then there exists a constant C2 = C^m, a, a±) such that 



0{t-°) + 0{ 1 k k- 2+l3 )- 



3) if a\ = a c , then there exists a constant C c = C c (m, a, o>%) such that, for any v € (0, 5)) 

D k {t) 



- C c u c {k) 



0{t-i +v ) 



(1.7) 



(1.8) 



(1.9) 



uniformly in k. 



Where u c (k) = I t k l e 1 ~ t dt and (3, 7, 9 and \i are given in fl.5\) . 
Jo 

Remark 1.2 The integral u c (k) — / t k ~ 1 e~ TTr[ dt can be rewritten as 

Jo 



u c {k) ■ 



k-2 k-2-; 

E E 



i=0 1=0 

'fe-1 



k - 1\ (k - i - I - 1)! 
i J (k - i)\ 



^ -^\k—i—l—l 



E 

i=0 



(k-i)\ 



+ 00 



r 2 e- pt dt. 



With help of computer calculation, u c (k) satisfies 

lim lnu c (fc)/(— k) = lim (— lnfc)/lnu c (fc) = 0. 



k — >c 



Based on Theorem ll.il we can obtain following two corollaries, which provide a complete distinction 
with respect to the parameters between the degree sequences for the present model. 
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Corollary 1.2 If the parameters satisfy that 

1) a > 2/3; or 

2) a < 2/3 and ot\ < a c , 

then the present random graph process has the power law degree sequence {1. Tjj . 



Corollary 1.3 Assume a < 2/3. 

1) If a c < ol\ < 2a c , then the present random graph process has the exponential degree sequence 



2) If a\ = a c , then the present random graph process has the critical degree sequence ll.9\) . 

Remark 1.3 When a > 2/3, for any ol\, the inequality ct\ < a < a c = Aa — 2 holds always, therefore, 
the part 1) of Corollaru \1.2\ follows from the part 1) of Theorem \l.l\ The part 2) of Corollary \1.2\ and 
Corollary \1.3\ are straightforward from Theorem 

Remark 1.4 A special case of the part 1) in Corollaru \1.2\ is a = a\ = 1. In this case, the model 
has a power law degree sequence as Ck~ 3 , which coincides with the result of [III. Furthermore, for 
any a € (1/2, 1] and ot\ = 2a — 1, the model has the degree sequence Ck~ 3 . 

Remark 1.5 The results are unclear for the following case: a < 2/3, 2a c < a.\ < a. Clearly, this 
case can only appear when a < 4/7. It is natural to conjecture that the model possesses an exponential 
degree sequence in this case. 

The methodology of the proof for the main results follows the standard procedure which can be 
found in [12] and [13] ■ The rest of the paper is organized as follows. In Section 2, we bound the degree 
of vertex in G t . In Section 3, we establish the recurrence for Dk(t) and then derive the approximation 
of Dk(t) by a recurrence with respect to fc. Finally, in section 4, we solve the recurrence in k using 
Laplace's method 18j and finish the proof of Theorem ll.il 




2 Bounding the Degree 

For times s and t with 1 < s < t, let d Xs (t) be the degree of vertex x s in G t . If x s is not added in 
Time-Step s, i.e., at Time-Step s, one of the other two substeps is executed, put d Xe (t) = 0. In this 
section, we will concentrate on the upper bound of d Xs (t). 

For the present model, the estimation for vt is derived in [13] as 




\v t — a\t\ < ci 1 / 2 logt, qs, 
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for any constant c > 0. We say an event happens quite surely (qs) if the probability of the compli- 
mentary set of the event is 0(t~ K ) for any K > 0. 

For the estimation of et, it can be derived by the same argument as in 13J that 

\et-rft\ <ct 1/2 logt, qs, (2.1) 

for any constant c > 0. 

By a standard argument on large deviation (see e.g. [21] and [23]), one further has: for any e > 0, 
there exists c\ , c% > such that 

F(e t < (r? - e)t) < ci exp{-c 2 t}, (2.2) 

for all t > 1. 

The following is our bounding for d Xs (t), note that our result is based on the exact estimation 
(|2.2[) for et. In our opinion, to bound the degree of vertex effectively, aforehand good estimations for 
et are necessary. 

Lemma 2.1 For any a G (1/2, 1] and a\ S (0, a], 

d x .(t) < (0O Pe (logi) 3 qs, (2.3) 

where p e is given in (|1.4|) . 

Proof: Fix s < t, suppose that x s is added in Time-Step s. Let X,- = d Xs (r) for r = s, s + 1, . . . , t 
and let 

ATMe(logi + l)' 1 ; 

12m 2 

where JV be large enough and will be determined later, and M t — . Let Y be the {1, 2, 3}-valued 

r\ — e 

random variable with P(Y = 1) = oti, P(Y = 2) = a — ot\ and P(Y = 3) = 1 — a. Then conditional 
on X T = x and e T > to, we have 

X XX 

X T+1 = x + I {Y= i)B{m, —) + I {Y = 2 }B{2m, —) - I {Y = 3} S(m, — ), (2.5) 

x 

where B(m,p) is the Binomial random variable with parameter (m,p) and S(m, — ) is the super 

e r 

geometric random variable with parameter (e r , x, m). 

Noticing that A is small enough for large N, using the basic inequality 



e y < 1 — y + 2y 2 for small y > 



G 



and the fact that S(m, — ) < to, (|2.5[) implies 



E (E(e AXT+1 | X T = x, e T ) | e T > to) 

+ (a - cti) 



< e Xx i ax 



2m 



-(1-a) 



1 - AE ( ^(m, — )J + 2A toE ^5(m, — ) 



Using the inequalities 



and 



e v < 1 + y + 2y 2 for small y > 



(1 + t/) m < 1 + my + —y 2 for small y > 



to the right hand side of (|2.6[) in turn, we get 



E (E(e A * T+1 | Jf T = x, e T ) | e r > to) < e Aa: <{ 1 + [(a c - ax) + 12mA 



mXx 



2e T 



< e Xx i 1 



< e 



mAx 
2e T 



max < (a c — ax) 



r\ — e 



12toA 



. I mAx max {(a c — ail, — -j , 

^ 1 ~t~ ^ m J (l + M e A) 



< exp < Ax 



i max{(a c -ai),^} 



2e T 



-(l + M e A) 



Now, we express E(e T+1 | X T — x) as 



E (e AXT+1 l T =i)=E (E(e AXT+1 | X T = x, e T )) 
= E (E(e A ^ +1 | X T = x, e T ) \ e T < to) P(e T < to) 

+E (E(e AXT+1 | X T = x, e T ) \ e T > to) P(e T > to) 
=:/ + !/. 

On one hand, conditional on X T = x and e T < to, X t+ x <x + m<e T + m< 2m holds always 

I < e 2m P(e T < to). 



On the other hand, II can be expressed as 

II = E (E(e AXT+1 \X T =x, e T ) \ e T > m, e T > (77 - e)r) 
xP (e T > to, e T > (77 — e)r) 

+E (E(e A ^ +1 I Af T = x, e T ) | e r > m, e r < (77 - e)r) 
xP (e T > to, e T < (77 — e)r) , 
by (|2.8p and the fact that x < e T , 

max{(a c - c*i), ^}to 



77 < exp < Ax 



2(r ? -e)r 



-(1 + M e A) 



+e X(,-e). exp max{(a e -g 1 ),^}mA (i + < _ 



< 



exp • 



I Ax l + ^(l + M e A) } + C"e A ^- £ ) r P(e r < (r? - e)r) 



(2.11) 



for some constant C" = C'(a,ai,e,m) > 0. By (|2.2|) and (|2.4p . choosing A~ large enough, then there 
exists constants C3 , C4 > such that 



77<expjAx 1 + — (1 + M e A) | + c 3 exp{-c 4 r}. 



(2.12) 



Combining (|2.9[) - (|2.12[) . using (|2.2p again for (|2.10|) , then there exists constants c 5 ,c 6 > such 



that 



Thus 



E( 



(e AX - +1 I X T = x) < exp |Ax 1 + ^(1 + M e A) j + c 5 cxp{-c 6 r}. 



E(e A * T+1 ) < E 



exp {* T A (l + PJ±±M^ }) + C5 cxp{ - C6T} . 



(2.13) 



(2.14) 



Now, put At = A and A r _i = A T (1 + Edl+Mi^zly Obviously if A s is small enough, then (|2.14|) 
holds for A r+1 , r = s, s + 1, . . . , t — 1. This will imply that 

t 

E(e AXt ) =E(e AtXt ) < E(e A * x =) + c 5 ^ exp{-c 6 r} < e mAs + C" (2.15) 

r— s 

for some constant C" > 0. 



Let A = 
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— , note that A can be taken small enough uniformly in t by taking A" large 
NM e {\ogt + 1) 

enough. Now provided A r < A, we can write 



A T — 1 < Ar 1 



p e (l + M e A) 



and then 

t 



Xs < A J|h + Ml±Mj<ioA(V S r« 

r—s ^ ' 

which is < A by the definition of A. 

Put u = (t/s)P'{\ogt) 3 , by (|2~T5j) we get 

P(X t >u)< (e mA * + C")e- Xu = 0{r K ) 

for any constant K > and the Lemma follows. □ 

Remark 2.1 For any n large enough, p e can be retaken as max < — ^— , — \, in fact, this can 



2(77 — e) ' n 
( 2(7/ - e) 1 

be done by enlarging a c — ot\ to max < (a c — ai), > instead in \2. 71). Thus, in the case of 

y nm J 

«i > o: C; p c can fee taken as 1/n. Certainly, if this is done as above, constants M e and N should be 
retaken correspondingly. 

3 The recurrence for D^t) 

In this Section, we follow the basic procedures in [13] to establish the recurrence for Dk(t). Put 
D-i(t) = for all t > 1. For k > 0, we have 

D k (t + l)=D k (t) 

+ (2a - oxJmE ( + VL^l® + 6t > ) P(e t > 0) 

Hi-«)M {k+1)Dk+l{t) - k -^+o(- 

\ e t e t V e * 



e t > m I P(et > to) 



+ai/ fc=m P(e t > 0) + 0(P(e t = 0)) + 0(P(e t < to)). (3.1) 

Here A t denotes the maximum degree in G t and the term O ( — -J accounts for the probability that 

V et ) 

we create larger than one degree changes for some vertices at Time-Step t + 1. By (|2.2p and Lemma 



2. 11 we have 



^ < O^'-^logt) 3 ), qs. (3.2) 
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The term E 



kD k (t) 



E 



f kD k {t) 



\ e t 

kD k {t) 



E 



et > Oj can be expressed as 
e t > 0^ 

e t - < i 1/2 logt ) P f|e 4 - r7*| < t 1/2 \ogt e t > 



+E ( ^ D ^ \ et - v t\ > t^ 2 logt,e t > o) F(\e t - rjt\ > t^ 2 \ogt \ e t > 0) 
E (fcgfcft) | |e t - yt\ < t 1 / 2 logt) P(|e t - yt\ < t 1 / 2 logt | e t > 0) 

x(l + O^- 1 / 2 log*)) + 0(P(|e t - rjt\ > i 1/2 logt | e t > 0)), 



(3.3) 



where we used the fact that kD k (t) < 2et to hand the second term. For k > 1, we have D k (t) 
E(D k (t) | et > 0)P(e t > 0), so 

E(kD k (t) | |e f -tjt\ <t 1/2 logt)P(|e t — <* 1/2 logt | e t > 0) 

= kD k (t)-E(kD k (t) \\e t -r)t\ > t 1/2 logt, e t > 0) 

xP(|et-^l>^ 1/2 logt|et>0) 

= fcD fc (t) + 0(i • P(|e t - 77*| > < 1/2 logt | et > 0)). 



Thus, using (|2.1[) , we have for fc > 



Similarly, 



E 



E 



< : i 



kD k (t) 



e/ 



„) = *M)+0(t-Vaiogt). 

70 



et >m)=^ + 0(t-^lo gi ). 



(3.4) 



(3.5) 



(3.6) 



Substituting (|3.2p . (|3.5|) and (|3 .6[) into (|3.1D . using (|2.2[) again to the other terms, we derive the 
following approximate recurrence for D k (t): D-\(t) — for all t > and for fc > 

D k (t + 1) = D k (t) + (A 2 (k + 1) + B 2 ) Dk+ ^ t] + (A 1 k + B 1 + l) Dk " ' 



t 



t 



+ (A (k - 1) + B ) Dk - l{t) + ai I k=m + O^-^logi) 3 ), 



(3.7) 



where 



a I - a 

A 2 = ; Ai 

2a -V 



2 — cei , 2a — a-i „ „ , „ 

= — tt; B 2 =B = and Si = -1. 



2(2a-l) 2(2a-l) 
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Note that the hidden constant, write as L, in term 0(t Pe 1 (logt) 3 ) of (|3.7[) is uniform in fc, which 

follows from the fact that e t = 0(t) and kD k (t) < 2e t = 0(t) uniformly in k. 

Dk(t) 

If we heuristically put d k = — - — and assume it is a constant, we get 

dk = {M{k + 1) + B 2 )d k+1 + (A t k + B 1 + 1)4 

+ {A (k - 1) + B )d k -i + ai/*=m + O^-^logt) 3 ). 

This leads to the consideration of the recurrence in fc: d-\ = and for k > — 1, 

(A 2 (k + 2) + B 2 )d k+2 + + 1) + Bi)dk+i + (A fe + Bo)d fc = -c*i4 =ro _i. (3.8) 

The following Lemma shows that, on certain conditions, (|3.8j) is a good approximation to (|3.7|1 . 
Note that our Lemma is a generalization of Lemma 5.1 in [13j . 

Lemma 3.1 Let dk be a solution for LS. 8\) such that \d k \ < ^ /or k > and a constant C . We have 
1) if a i < a c , then, for any v G (0, 1 — p c ), i/iere exists a constant Mi > smc/i that 

\D k (t)-td k \ <M l t"'+\ (3.9) 

/or all t > 1 and fc > — 1; 

i/a c < ai < 2a c , iften there exists a constant M 2 > such that 

\D k (t) - td k \ < Aht 1 - 9 , (3.10) 

for all t > 1 and fc > — 1, where 9 is given in \1.5\) . 

Froo/. Let fc (i) = D k (t) - td k and fc = k a {t) = \t p ' (log tf\ . Lemma O implies 

< I?fc(t) < i~ 10 for k > k [t). (3.11) 

Proof of part 1): Equation 13 . 1 1 [> and d k < C/k imply that (|3.9j) holds for fc > fco uniformly, i.e., 
there exists a constant Ni > 0, independent to fc and r, such that 

\D k {t)-td k \ = \e k (t)\<N 1 t"' 

for all fc > fco(^) and t > 1. 

Recall that the hidden constant in 0(i Pe_1 (logr) 3 ) of ()3.7p is denoted by L. For any v e (0, 1 — p e ), 
let i? > L satisfying 

Lt^'- 1 (log*) 3 < Rt p ^ +U - 1 
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for all t > 1. Let N 2 = — ^- + 1, take a > such that 

l-^r-(l + <7)(l-p e -f)>0, 

and take <5 G (0, 1) such that 

S 1+a < e- 1 < S. 



Let t\ > be an integer such that 



, , . -1 2(2a-l) 
/Li 2 — ai 



^<d i v +i , fi-i^a™<* 



and 

^ 1+CT < f 1 - , - , 

t+lj \ t+l 

for all t > h and I > N 2 . 

Now, for the above t x , let > N x satisfying 

|©fc(*)| < N z t p ' +y for all 1 < t < ti and k > -1. 

Take 

Mi = max{N 2 ,N 3 }. 

We will prove that (|3.9[) holds for the above Mi by induction. Our inductive hypothesis 

Hi ■ |@fc(*)| < M 1 t p - +U for all k > -1. 

Note that (pTfl)]) and ([XTT]) imply that ft| holds for 1 < t < t x . 
It follows from and flJTHD that 

e fc (t + 1) = e*(*) + A 2 (k + i)^±iW + {Al k + b 1 + i)^ 

+ A o(fc - + (^-i(logt) 3 ). 

For i > ij, by (pHI)) . we have t + A\k + B x + 1 > and then (pTTg)) implies 
|9*(t + 1)1 < A 2 (fc + l) |9fc + lWI + (i + Aifc + Si + 

< (t + ^ 2 (fc + 1) + A x k + B l + 1 + A (k - l))M x t p * +v - 1 + Rt»^ +V 
= (t + A 2 + B 1 + l- A )M 1 t p ' + "- 1 + Rt p < + »-\ 
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Let Eq = A 2 + Bi + 1 — A ~ (aii — a c )/a c , noticing that ai < a c , we have Eq < 0. Then, combining 
1(532)1 . dHUH) and (J3H7J), we have 

(i + eo^/ii^- 1 + RtP^- 1 M x tP^ v + RtP- +u - 1 
Mi(t+l)P'+ v ~ Mi(t+\)P'+ V 

l-R/Mj _(i_ pe _„) 
t + 1 f / 1 \ t+1 1 * + l 



< 1. 

The induction hypothesis has been verified and the proof of part 1) is completed. 

Proof of part 2): In this case, we have a c < cei < 2a c and then, for some v s (0, 1/2), £o < p e + v < 
1 — 9 ( note that in this case p e = 1/2). Same as what we have done for part 1), for certain a > and 
S 6 (e -1 , 1), we have 

for sufficient large t and M 2 . This is enough for a inductive proof of (|3 . 1 0[) . □ 

Remark 3.1 [R emar k 5.2 in [13]] Lemma\3J\ implies that if there is a solution for \3.8\) such that 
dk < C/k, then lim^oo D k (t)/t exists and equals to d k - In particular, it is shown that: if there exists 
a solution for \3. 8]) such that dk < C/k, then the solution is unique. 

4 Solving (13.81) and the proof of Theorem 11.11 

In order to solve (13. 8p . let us consider the following homogeneous equation 

(A 2 (k + 2)+B 2 )f k+2 + (Ai(k + l)+Bi)f k+ i + (A k + B )f k = 0, k>l (4.1) 

which is solved by Laplace's method as explained in (18) . 

For k > 1, we construct function f k has the following form 

,b 

fk = / t k ~ x v(t)dt, (4.2) 



where constants a and b, and function v(t) are to be determined later. 
Integrating by parts 

,6 

kf k = [t k v(t)] b a - / t k v'(t)dt. (4.3) 
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Let 

4>i(t) = Mt 2 + Ait + A , 0o (t) = B 2 t 2 + Bit + B 
Substituting (|4~2j) and (|4~3]) into (|4~T|) . we obtain 



6 

fe-i. 



[t k Mt)v(W a - t k Mt)v'(t)dt+ I t k - 1 <t> Q {t)v{t)dt = Q. (4.4) 

J a 

Equation (|4.1j) will be satisfied if we have 



V '(t) Mt) 



(4.5) 



v(t) t<j) X {ty 

and 

= 0. (4.6) 
Let a = and & equal to a root of v(t)4>±(t) = 0, the parameters a and b can be determined satisfying 

Obviously, 0o(^) an d 4>i{t) can be rewritten as 

Mt) = -t; 4>i{t)=At 2 -{A + B)t + B = A(t-l){t-B/A), (4.7) 

where 

A = i-4, B = 2a - ai =A+ ac ~ ai . (4.8) 
2a — 1 a c a c 

Now, we solve the equation (|4.ip in the following cases: 1), ai < a c ; 2), ai > a c and 3), «i = a c 

respectively. 

For case ai < a c , we have B > A, then the differential equation (|4.5p is homogeneous and can be 
integrated to derive 

v(t) = (t-lf(t-B/A)-P, (4.9) 

where = 1/(B — A) is given by (|1.5[) . 

Since in this case /3 > 1, so by (|4.7p . the equation 

w(i)0i(t) = A(t - l) 1+/3 (t - B/A) 1 -' 5 = (4.10) 



has a unique root 1. Thus, the parameter 6=1 satisfies (I4.6p . 

Substituting the parameter b and the function v(t) into (14. 2p and removing a constant multiplicative 
factor, we obtain a solution ui(fc) to (|4.ip for fc > 1: 

ui(k)=j\ k - i (j^ydt, (4.ii) 

where £ = A/5. 

The order of the function ui(k) with respect to fc is given by the following Lemma. 
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Lemma 4.1 [Lemma 6.1 in [13]] Let k > L rften 



ui(fc) - (1 + 0(fc- 1 ))X'i/c- (1+ ' 3) (4.12) 

/or X>i = T>i(a, a±) a fixed constant. 

In Case of a± > a c , we have B < A, and equation (14. 5|) has the same solution as (|4.9p . In 
addition, under the conditions (| 1 . 2|) and (|1.6|) . one further has /? < — 1, and then the equation (|4.10|) 
has a unique root 7 := as given in (|1.5p . So we can take 6 = 7 to satisfy (|4.6p . Thus 

" (t >-r^ i (^r*-^r , "(££r* 

is a solution to (|4.ip for A: > 1. 
By Lemma 14. 11 we have 

u 2 (fc) = (1 + Oik-^V^k- 1 ^ (4.13) 

for some fixed constant T>2 = ^2(0;, a±). 

Finally, we consider the case of cti = a c . In this case A = B and the equation (|4.5p can be 
integrated to derive 

v (t) = e"^ 1 -') 

with fj, — l/A given in (jl.5p . With same argument as in cases 1) and 2), take b = 1 and define 

u c (k) = [ t^e-^dt, 
Jo 

then u c is a solution to (|4.ip for /c > 1. 
Crudely, 

u c (k) < [ t k - x dt = l/k. (4.14) 
Jo 

The precious representation of u c (k) can be found in Remark 1 1.21 

Note that in all the three cases, u\, U2 and u c do not satisfy equation (|4.ip when k = 0. In fact, 
as calculated in [TSJ, for i = 1,2 or c, we always have 

2A 2Ui {2) + (A, + BiH(l) = [0i(t)«(t)]J; - -0i(O)«(O) + 0. (4.15) 

Now, we are going to solve p.8p . By Remark 13.11 we only need to construct a solution for (|3.8p 
which satisfies the requirements of Lemma 13.11 Actually, we will construct such a solution based on 
the solution of (|4.ip given above. 

Denote by g the solution for (|4.ip . i.e., g = U\,U2 or u c in the three cases respectively. 
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For m > 1, define Wk = for fc > m, u> m -i = — ai/[(m — 1)A)] an d for j = m — 2, m — 3, . . . , 1, 
let Wj be such that 

A 2 {j + 2)w j+2 + {Ai{j + 1) + B{)w 0+1 + A jwj = 0. 

Then Wk satisfies (|3.8[) for fc > 1. Therefore, any linear combination of g and w is a solution of (|3.8|) 
for fc > 1. 
Now, let 

D= 2A 2 w 2 + (^! +gi)wi ^ A 2 (£>g(l)+wi) 
2A 2ff (2) + (Ai + Bi)s(l) an Si 
Note that S and d depend on g = ui, u 2 and n c respectively. By (I4.15[) . D is well-defined. 

Define 

'0, if k = -1 

<4 = < d, if fc = • 

Dg(k)+Wk, otherwise 

It is straightforward to check that dk given above is the solution of (|3.8p , by (|4.12p , (|4. 13|) and (|4. 14(1 , 

we know that dk satisfies the requirements of Lemma 13.11 

For m = 1, we can take 

n = 01 SAgffCl) 

2A 2 « 7 (2) + (A 1 +S 1 ) 5 (1)' S x 

and directly define 

!0, if fc = -1 

d, if fc = 

Dg{k), otherwise 

Similarly, in this case dfc is also a solution to (|3.8| which satisfies the condition of Lemma 13.11 

Proof of Theorem ] By the construction of the solution dk and Lemma 13.11 the theorem follows 
immediately by taking 

J2A 2W ^(A 1+ B 1)W1 )V iQrm>i 

2A 2Ul {2) + (Ai + fli)ui(l)' 
d={ fori = l,2,c, 

«iA f = 1 

2A 2 u 4 (2) + (Ai+B 1 )u i (l)' 01 ™ 

where Pi and £> 2 are given in P~T2")) and (07T3J), P c = 1. □ 
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